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We calculate the transport properties of three-dimensional Weyl fermions in a disordered environ¬ 
ment. The resulting conductivity depends only on the Fermi energy and the scattering rate. First 
we study the conductivity at the spectral node for a fixed scattering rate and obtain a continuous 
transition from an insulator at weak disorder to a metal at stronger disorder. Within the self- 
consistent Born approximation the scattering rate depends on the Fermi energy. Then it is crucial 
that the limits of the conductivity for a vanishing Fermi energy and a vanishing scattering rate do 
not commute. As a result, there is also metallic behavior in the phase with vanishing scattering 
rate and only a quantum critical point remains as an insulating phase. The latter turns out to be 
a critical fixed point in terms of a renormalization-group flow. 

PACS numbers: 05.60.Gg, 66.30.Fq, 05.40.-a 


I. INTRODUCTION 

Since the discovery of the fascinatingly robust transport properties of graphene [IHi , there has been an 
increasing interest in other two-dimensional systems with similar spectral properties, such as the surface 
of topological insulators In all these systems the transport is dominated by a band structure, in 

which two bands touch each other at nodes. If the Fermi energy is exactly at or close to these nodes, 
the point-like Fermi surface and interband scattering lead to particular transport properties, such as 
a robust minimal conductivity. Based on these results, an extension of the nodal spectral structure to 
three-dimensional (3D) systems is of interest [iM3- In 3D the Fermi surface is a sphere with radius \Ep\ 
rather than the circular Fermi surface in 2D, which is either occupied by electrons {Ep > 0) or by holes 
[Ep < 0). For Ep = 0 the conductivity vanishes in the absence of impurity scattering in contrast to the 
minimal conductivity of the 2D system. On the other hand, sufficiently strong impurity scattering leads 
to a conductivity at the node Ep = 0. Thus, an important difference between 2D and 3D Weyl fermions 
is that there exists a metal-insulator transition in the latter, which is driven by increasing disorder 
[ttI [29l - [3^ . This transition is similar to the metal-insulator transition caused by decreasing random gap 
fluctuations in a system of 2D Dirac fermions (^. Is^. On the other hand, it is quite different from an 
Anderson transition from a metallic state at weak scattering to an insulating state at strong scattering, 
which is caused by Anderson localization at strong scattering [s^, [111 . 

There is agreement between the various approaches, based on self-consistent, perturbative and nu¬ 
merical methods, on the existence of a transition from a 3D Weyl semimetal at weak scattering to a 
diffusive metallic behavior at stronger random scattering [lol - [^ . This transition can be characterized 
by a vanishing density of states at the Weyl node (i.e., the scattering rate or the imaginary part of the 
self-energy) and a nonzero density of states in the diffusive phase. However, the transport properties for 
the weak scattering regime are still under discussion. In particular, a recent study indicates that there 
is a metal-metal transition rather than a insulator-metal transition for 3D Weyl fermions with a critical 
point [sll . We will address this problem in the subsequent calculation, using a weak scattering approach 
(WSA). 

Calculations of quantum transport consist usually of two steps: Determining the scattering time (or 
scattering rate) within a self-consistent solution of the Dyson equation, also known as the self-consistent 
Born approximation (SCBA), and determining the conductivity by a self-consistent solution of the Bethe- 
Salpeter equation (BSE). This approach, in particular the solution of the BSE, is rather complex due 
to the existence of many modes. Not all of them are relevant for the transport properties because some 
decay quickly. From this point of view it is easier to project at the beginning only onto those modes which 
do not decay quickly but control the transport properties on large scales. For a system with spectral 
nodes in a disordered environment these modes are a result of a spontaneously broken chiral symmetry 
[sgl - ld^ . We will employ this idea here to 3D Weyl fermions in order to calculate the conductivity. For 
this purpose it is important to identify the underlying symmetries of the two-particle Green’s function. 
Then spontaneous symmetry breaking is characterized by its non-vanishing order parameter which is the 
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scattering rate. There is a metallic phase with long range correlations (i.e. diffusion), whereas in the 
insulating phase the symmetry remains unbroken. 

The paper is organized as follows: In Sec. II we define the model and discuss the symmetry properties 
of the two-particle Green’s function. Then the DC conductivity is calculated within a weak scattering 
approach (Sect. Ill), using an expansion in powers of the scattering rate rj. This provides us a formula 
for the DC conductivity, which is discussed in Sect. IV at the node (Sect. IVa) and away from the node 
(Sect. IVb). Our discussion includes a comparison with the results of the Boltzmann approach and with 
results from an approach based on the SCBA and BSE of Refs. [H, 


II. MODEL 

The three-dimensional Weyl Hamiltonian for electrons with momentum p is expanded in terms of Pauli 
matrices r, {j = 0,1, 2, 3) as 


H = vft-p—Uto with t = {ti,T 2 ,T 3 ) 


( 1 ) 


with Fermi velocity vp- U is a disorder term, represented by a random potential with mean (U) = Ep 
(Fermi energy) and variance g. The average Hamiltonian {H) generates a spherical Fermi surface with 
radius \Ep\, and with electrons (holes) for Ep > 0 {Ep < 0). Physical quantities are expressed in such 
units that vpfi = 1. 

The electronic conductivity, obtained as the response to a weak external field with frequency a; ~ 0 

'141141 


c^(w) = X! > (2) 

r 

is given by the correlation function of the Green’s functions {H ± z)“^ 

Arr'{uj) = iirR{Tr 2 [{E[ + z)~^,{E[ - z)~,l]) with z = uj/2 + ie, (3) 

where (...) represents disorder average and Tr 2 is the trace with respect to the Pauli matrix structure. 
This expression, often called the two-particle Green’s function, has two different energies ±z for the same 
Hamiltonian H to create two independent Green’s functions {E[Ez)~^. Now we represent this two-particle 
Green’s function by two different Hamiltonians and one energy: We define the pair of Hamiltonians 

H± = Pin + P2T2 ± {psn -Utq) , ( 4 ) 

where Hp = H. The matrix transposition ^ relates Hp and i4_ through the identity 

TiHln = -Hp , (5) 

since pj = —pj. This allows us to write for the correlation function ([3]) 

Arr'{uj) = - lim(Tr2[(I4+ -f z)~pTi{H'^ + z)~pTi]) . (6) 

£—>•0 

Instead of two different energies Ez and the same Hamiltonian El, the two-particle Green’s function 
has now the same energy z but different Hamiltonians, namely Hp and H^. The relation ([5]) and the 
representation ([6|) reveals an internal structure of the model which leads to the Hamiltonian 


The Green’s functions {Hp + z) ^ 
element of the Green’s function {H -b z)~^ 



(Hp 

0 

0 


0 

iJ_ 

0 

H = 

0 

0 



\ 0 

0 

0 

and {H^ + 

z)-i 

in Eq. 


0 

0 

0 


(7) 


are just the first and the third diagonal 
This indicates that the transport properties of the original 
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Hamiltonian H, which requires two different energies ± 2 , are related to the transport properties of the 
extended Hamiltonian H at the same energy z. 

The extended Hamiltonian H, its symmetries and its relation to diffusive transport were studied pre¬ 
viously [H, 113, m, . In particular, it was found, together with property (O, that the matrix 


/ 0 0 (pin 0 \ 

0 0 0 (P 2 T 1 I 

p'lTl 0 0 0 I 

V 0 P2'Ti 0 0 / 


( 8 ) 


with scalar variables cpp (/?'- anticommutes with H: SH = —HS. This relation implies a non-Abelian 
chiral symmetry [^, |4^ : 


e^He^ = H 


( 9 ) 


which is a symmetry relation for the extended Hamiltonian H in Eq. ©• The term proportional to z 
in the Green’s function G'(z) breaks this symmetry due to ^ 1 , and therefore, lim 2 ^o[G'(z) — G'(—z)] 
plays the role of an order parameter for spontaneous symmetry breaking: 


G(z) - G(-z) = - 2 zG(z)G(-z) = - 2 z{H'^ - z^)-^ . 


( 10 ) 


Since the diagonal elements of this expression are proportional to the density of states at the node when 
we take the limit z —>■ 0, a non-vanishing density of states indicates spontaneous symmetry breaking. 
The role of a non-vanishing density of states at the node as an order parameter for a diffusive metallic 
phase was also discussed in Refs. [H, [H, IMj. 

Following the recipe of Ref. the correlation function ([S]) can be expressed as a diffusion propagator. 
This is used in the next section, where we focus on the long-range behavior of Arr>(uj) to calculate the 
conductivity < 7 ( 0 ;). 


III. WEAK-SCATTERING APPROACH 

The scattering rate 77 is defined by 

^ 2 Tr[{{H + ie)-^){{H-ie)-^)] ' ^ ^ 

This definition is motivated by the assumption of a complex self-energy for the average one-particle 
Green’s function {{H+ie)~^) , whose imaginary part is the scattering rate (cf. App. |^. The corresponding 
scattering time r is r = h/rj. r] can either be calculated, for instance, within the SCBA I, m 113 or 
it can be taken from experimental measurements. As discussed in the previous section, a non-vanishing 
scattering rate indicates spontaneous symmetry breaking. Since the broken symmetry is continuous, 
there exists a massless mode. The latter is reflected by the relation 

V Arr'ioj) = Tr2{[{H + uj/ 2 )-\H - u:/ 2 )-\r) = --Tr2[{{H + uj/ 2 )-^) - {{H - w/2)-i)] (12) 

UJ 

r 

which diverges for a vanishing symmetry breaking term w ^ 0 due to long-range correlations. For the 
correlation function a similar but more elaborate analysis yields a diffusion propagator (sJl, whose 
Fourier components read 


Ag{uj) 


V 1 
g iujl2 + Dq^ 


(13) 


This agrees with (HU for g = 0 when we use the self-consistent approximation {[H + ie + zo/2)Ki —Y^jg 
of App. 1^ This is not an accident but a consequence of the fact that the self-consistent approximation 
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represents the saddle point of the corresponding functional integral [d^. The prefactor of the term is 
the diffusion coefficient 

n 91 f rr. (d{{H) + ir])-^ d{{H) - iiq) 

= - SkT 



Thus, the DC limit w ^ 0 of the conductivity formula in (I5|) and the correlation function (IT5)l reproduce 
the Einstein relation 


which gives with the right-hand side of Eq. (d for 3D Weyl fermions the integral 

2 il'^ + k'^)^ + Ej,{2ri‘^+ 2k^/3 +Ej,) k'^dk 

[(^ 2_£;2 + fc 2 ) 2 + 4 ^ 2 ^ 2]2 27 r 2 


(15) 


(16) 


with momentum cut-off A. Thus, the conductivity depends on the disorder strength g only through the 
scattering rate rj. 

A diffusion propagator can also be calculated from the BSE, as demonstrated recently for 3D Weyl 
fermions [^,|^. However, the derivation of the propagator (fT^ from the symmetry ([5]) has the advantage 
that it is simpler and that we obtain the diffusion coefficient D in (fT4l) directly as a quadratic form of 
Green’s functions. 


IV. RESULTS AND DISCUSSION 

In the following we present and discuss the results which are obtained from the conductivity (T(ry, Ep) 
in Eq. (dll). This expression is subtle in the limit of a vanishing scattering rate 77 , since the latter appears 
as r]^ in front of an integral that diverges for 77 —>■ 0. This makes the conductivity very sensitive to the 
order of the limits Ep —?> 0 and 77 —>■ 0 in the case when the scattering rate vanishes at the node. Since the 
conductivity depends on 77 and Ep separately, we consider first properties exactly at the node Ep = 0, 
where results are simple, and then the more complex results when 77 depends on Ep. For the second part 
we employ the SCBA to determine the function r]{Ep) and calculate the corresponding conductivity. 


A. Transport at the spectral node 


At the node {Ep = 0) the DC conductivity in Eq. (IT6)) is reduced to the expression 


2 k'^ ^k 

Jo (772 + P)2 ^ 


9 

e T] 
2Tr^h 


arctan(l/^) 


c 

1 + C^ 


(C = ??/A) , 


which becomes for 


a ~ 


e 


2 


Airh 


(17) 


(18) 


In contrast to the 2D case, where a = l-Kh, the 3D case gives a linearly increasing behavior with 
respect to the scattering rate. This result was derived directly (i.e., without using Eq. (fTCl) ') by Fradkin 
some time ago M- With a disorder dependent scattering rate he also obtained a transition for a critical 
disorder strength g^, where the conductivity vanishes for g < gc and increases linearly for g > gc- 

The linear behavior indicates an unconventional transport because in the classical Boltzmann approach 
for one-band metals the conductivity decreases with increasing scattering rate: ap = ne'^/mhg {n: 
electron density, m: electron mass) (46l| . This remains true when we include the band structure of the 
Weyl fermions in the Boltzmann approach: ap is nonzero at the node for any scattering rate and even 
diverges with vanishing disorder as 


1 e^Vph 
2tt g 


crp = 


(19) 
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where g is related to the density of impurities and the impurity potential Uq by the equation g = riiU^ 
[33 ■ The disagreement between the expressions in (fT71) and (fTOl) can be explained by interband scattering, 
caused by particle-hole creation processes, which has been ignored in the Boltzmann approach. On the 
other hand, the increasing behavior of (ED for small 77 turns into a decreasing behavior for larger 77 , as 
one can see in Fig. [I] indicating a crossover from quantum transport for weak scattering to conventional 
Boltzmann transport for stronger scattering. 

In this context it is also interesting to study the finite-size effects of the conductivity by considering a 
cubic system of finite length L. The /3-function /3 = dhia/d\ii.L describes the finite-size scaling of the 
conductivity. It can be calculated from Eq. dT51) by replacing the lower integration boundary with I/L, 
which assumes non-periodic boundary conditions. Then a{r], L) in Eq. (1171) is a monotonically increasing 
function for increasing size L with 

for Lt] ^ 1, where a* is the expression ED- The corresponding /3-function reads 

/3^2[l-a(77,L)/a*] , (21) 


which vanishes at the 77 -dependent fixed point a*. This differs from the 2D case only by different fixed 
points a*, where in 2D it is a universal constant a* = j'Kh (T7j | and in 3D it is the 77 -dependent expression 

ED- 

It should be noticed that a* is not a critical point because it is an attractive fixed point. But since for 
3D Weyl fermions cr* depends on the scattering rate 77 , we have a line of fixed points for 77 > 0. Thus, 
the endpoint tr* = 0 for 77 = 0 has the feature of a critical point because any change of 77 drives us away 
from this endpoint, as illustrated in Fig. [51 It indicates a transition from an insulator (77 = 0) to a metal 
(77 > 0). The transition is driven by increasing disorder, since the scattering rate is a monotonic function 
of the disorder strength g. rj is also the order parameter for spontaneous symmetry breaking ED, which 
can be calculated from Eq. ED within SCBA. From the solution of the self-consistent equation dMD we 
get for 7 = < 7 / 271 ^ and rj ^ 0 the linear behavior 


rj -(Ay - 1) ©(Ay - 1) (y = g/2Tr^) (22) 

TT 

with the step function 0. For y < yc = 1/A we have no spontaneous symmetry breaking. Thus, 77 as 
well as the DC conductivity vanish strictly. When we approach yc from above there is linear behavior for 
the scattering rate, which agrees with the numerical calculation of Kobayashi et al. [ 2 ^ . At the critical 
point itself we obtain from the Einstein relation (1151) a finite diffusion coefficient 


D{9c) 


9ce^ 

47r 


The results of the DC conductivity from previous self-consistent studies, based on a combination of 
SCBA and a self-consistent solution of the BSE [ 2 ^, [s^ , are summarized and compared with our results 
of the WSA in Table [H For sufficiently large scattering rates the Boltzmann approach, the solution of the 
BSE and the result of the WSA agree reasonably well, reflecting a rather conventional transport. This 
indicates that quantum effects, such as particle-hole pair creation, are dominated by impurity scattering. 
On the other hand, for smaller values of the scattering rate the conductivity exhibits a larger variety of 
results: The Boltzmann conductivity has a simple 1/g behavior, which is also found with the solution of 
the BSE in Ref. [ 2 ^, with a different constant prefactor though. In contrast, the approximative analytic 
solution of the BSE in Ref. has a characteristic dip down to zero at a critical gc and increases for 
g> gc and for g < g^: 


cr = (Til 1/g - l/5c|, (Ti = cr 


1 for g < gc 
3 for g > gc 


(23) 


When we compare this result with the WSA conductivity in Eq. (fTTl) it should be noticed that the latter 
was obtained by sending Ep ^ 0 first and then 77 0. As mentioned at the beginning of this Section 

the value of the conductivity depends on the way we take these two limits. Although nothing has been 
said in Ref. [s^ about the order of the two limits to get (l23l) . we will study in Sect. lIVBI the case when 
Ep and 77 go to zero simultaneously in Eq. (ED- Then we obtain a result similar to (1231) . 
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FIG. 1: The conductivity of 3D Weyl fermions Cl as a function of the scattering rate 77 at Fermi energy Ef = 0. 
It describes the crossover from the quantum regime for rj < 0.6 to the classical Boltzmann regime for 77 > 0.6. 



FIG. 2: The /3-function for different scattering rates, where the arrows indicate the flow toward the fixed points. 
The /3-function increases with 77 and creates a line of fixed points a* ~ rjc^/Ahn. This plot is based on an 
approximation near the fixed points according to Eq. 1211). 


B. Transport away from the spectral node 

The conductivity as a function of the Fermi energy is plotted at fixed scattering rates in Fig. [31 As 
we increase the scattering rate the effect of the node is washed out and the conductivity becomes flatter. 
This is similar to the behavior in Fig. [T] In other words, impurity scattering supports transport near the 
node whereas it suppresses it further away. Thus, we can distinguish a regime close to the node, where 
the conductivity increases with the scattering rate, and a more conventional regime further away from 
the node, where the conductivity decreases with the scattering rate, as also described by the Boltzmann 
approach. 


at the node 

Boltzmann approach [30] 

SGBA & BSE [2^ 

SGBA & BSE [3^ 

SGBA & WSA 

scattering rate rj 

0 

0 

2 (ff/5c - l)/7r 

‘2{g/gc - l)/7r 

conductivity a 

e^v%h/2ng 

Ae^VpIgh 

(Tie^/dvr/i 

Eq. (flTll. Fig. [1] 

away from the node g < gc 





scattering rate rj 

oc Ep 

hgEp/S-KV^ 

Epgc/{gc - g) 

Eq. (dll) 

conductivity a 

aB{0){l + 6El/E^) 

Ae^VpIgh 

(xEpgc/igc-g) 

Fig. a 


TABLE I: The scattering rate 77 and the conductivity at the node calculated with three different methods. The 
translation from Ref. is <? = ruu^ and the SGBA coefficient ci is given in Eq. (I23II . 
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Fermi energy 


FIG. 3: The conductivity (1161) as a function of the Fermi energy Ef at fixed scattering rates r] — 0.02 (full 
curve), rj = 0.05 (dashed curve) and 77 = 0.1 (dotted curve). 


So far we have considered the case that rj and Ep are independent. However, in general the scattering 
rate depends on 7 and Ep. For instance, the self-consistent calculation in App. based on the SCBA 
or saddle-point approximation, creates a scattering rate in Eq. (IA4I1 that depends on the Fermi energy: 


( 7 A -1) , 

/(7A-1)2 

2iEp 


I 7^77^ 

77r 


(24) 


The behavior of the conductivity in dT51) is affected by this result, since the limits ry —>• 0 and Ep —>■ 0 are 
not independent anymore. (1241) has two typical regimes, namely rj oc Ep near the critical point 7 A = 1, 
which leads to cr oc Ep, and rj oc Ep for 7 A <C 1, which leads to a non-vanishing conductivity for Ep —F 0. 
This implies that the conductivity vanishes only for 7 A ~ 1, whereas it nonzero above and below 7 A ~ 1. 
Thus, there is no insulating phase but only an insulating point for 7 A = 1, in agreement with the analytic 
result of Ref. [S^ . The behavior of the conductivity as a function of 7 A — 1 is plotted for different values 
of the Fermi energy in Fig. |4l It should be noticed, though, that the transport behavior is determined 
by the Ep dependence of the scattering rate of the special form in Eq. (l24l) . Using another form of the 
scattering rate as a function of the Fermi energy can lead to a substantially different behavior of the 
conductivity near the node. An example was observed in Ref. [30l| within a numerical solution of the 
SBCA and the BSE, where the scattering rate is exponentially small for 7 < 7 c. In this case a vanishing 
conductivity was found at the node also for 7 < 7c. 

V. CONCLUSIONS 

We have studied the DC conductivity of 3D Weyl fermions in the presence of random scattering. The 
relevant parameters in the conductivity (I16p are the scattering rate rj, which is an order parameter for 
spontaneous chiral symmetry breaking, and the Fermi energy Ep. Exactly at the node Ep = 0 there is a 
metal-insulator transition with a diffusive metal for 77 > 0 and an insulator for 77 = 0. The conductivity is 
linearly increasing with 77 up to a maximal value and decreases for stronger scattering rates, as illustrated 
in Fig. [TJ This non-monotonic behavior is in contrast to the constant conductivity in the corresponding 
2D system. It reflects the fact that the increased phase space of the 3D Weyl fermions suppresses the 
conductivity for weak scattering but also that stronger scattering implies a screening of the node such 
that the Boltzmann approach eventually becomes applicable. Further away from the node the behavior 
depicted in Fig. [3] agrees qualitatively with that of the 2D system which was also obtained in a 
quasiclassical approach [d^. The latter diverges as one approaches the node, which indicates that the 
full quantum approach is necessary near the node. 
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FIG. 4: The conductivity (I16II as a function of the disorder parameter A = 7 A —1 with a self-consistent scattering 
rate ? 7 (A) for different values the Fermi energy Ep- The curves tend to reach the critical point A = 0 by a sharp 
cusp as the Fermi energy approaches the node Ep = 0. Further away from the critical point the conductivity is 
much less sensitive to a change of Ep. 


The critical behavior at the node describes an unconventional phase transition for 7 = 7c which is 
driven by quantum fluctuations: In contrast to a conventional transition the symmetry broken phase 
with 77 > 0 is characterized by robust diffusion whereas the phase with unbroken symmetry (77 = 0 ) has 
a subtle behavior in terms of the conductivity because it is very sensitive to the limit Ep —^ 0. Thus, it 
is possible that we either have an insulating phase with vanishing conductivity when the scattering rate 
vanishes slowly with Ep or a. metallic phase when the scattering rate vanishes sufficiently fast with Ep. 
In the case of an SCBA calculation for 77 there is only a quantum critical point in the transport properties 
and metallic behavior above and below this critical point. It cannot be ruled out, though, that a different 
calculation of 77 leads to a different behavior. Thus, our discussion of the delicate limits 77 , Ep —>• 0 clarifies 
some of the contradicting results in the literature about the presence of a metal-insulator transition for 
3D Weyl fermions [^. Moreover, the fact that the conductivity depends only on 77 and Ep allows 
us to determine 77 independently with other approximations than the SCBA, and to insert the result into 
the conductivity (ITOl) . A possible step in this direction is a correction to the SCBA or perturbative 
renormalization-group approach in d — 2-expansion. The latter gives 77 (A) oc E^'^ [29l 1^. This result 
would lead to a vanishing conductivity below the critical point. 

In the regime 0 < 7 < 7 c disorder may also affect physical properties of Weyl fermions in another 
way. The reason is the existence of non-uniform solutions of the SCBA with an exponentially small 
contributions to 77 , similar to Lifshitz tails in the density of states of disordered systems [i^. In the 
case of 2D Weyl fermions this has been discussed in Ref. [s^. The problem, however, is always that 
the self-consistent equation is nonlinear and has many non-uniform solutions. Under certain plausible 
assumptions for the solutions 77 , the existence of exponentially small contributions for 7 < 7c has been 
discussed for 3D Weyl fermions in Ref. [^. The corresponding states might be localized then with no 
contribution to the conductivity at T = 0. However, for T > 0 thermally activated electrons may hop 
between patches of localized states and provide a hopping conductivity [^. Whether or not resonant 
tunneling without spontaneous symmetry breaking can occur in this case is an open question. 

Acknowledgement: I am grateful to David Schmeltzer for an extended discussion of Weyl fermions. 


Appendix A: Self-consistent approximation 

The first step is to study spontaneous symmetry breaking of the symmetry Q by a non-zero scattering 
rate within SCBA, following a similar approach as given in Ref. [^. The average one-particle Green’s 
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function then reads 


((i? + z)-^)«((i?) + z + E) 


-1 


(Al) 


and the self-energy is given by 


E = -g((i7) + E + ze)- 


(A2) 


The imaginary part of the self-energy E is a scattering rate rj. Then the self-energy reads in our case 
with the momentum cut-off A 


E = 70! 


A-flog 


o -I- A 
a — A 


For small Ep near the node we expand Eq. 

*(7-^-1) 


(7 = g/ 27 r^, a = Ep + E) . 
in powers of o up to second order to obtain 


E ^ —Ep + 


'(7A-1)2 2iEp 






The real part of E provides a shift of the Fermi energy: 


-E 


F + ^\ -- 


2 iEp 


77r 


= —Ep + e 


i-K j A. 


77r 


l 2 Ep 

77r 


where the sign is chosen such that we have a positive scattering rate. 

At the node Ep = 0 the self-consistent eq. (|A 3 p reduces to rj = rjI with 


7 = 7 [A — 77arctan(A/77)] . 

There are two solutions, namely rj = 0 and 7 7^ 0 with 

1 


A7 = 


1 — C arctan(l/<;^) ’ 


C = 7/A . 


(A3) 


(A4) 


(AS) 


(A6) 


A nonzero 77 reflects spontaneous symmetry breaking with respect to ([S]). Such a solution exists for 
only at sufficiently large 7. Moreover, 7 vanishes continuously as we reduce 7. A nonzero rj is proportional 
to the density of states at the Fermi level. However, even for 77 = 0 there can be a nonzero local density 
of states due to localized energy levels, which are not counted in 77 within the SCBA. For ^ ~ 0 we obtain 
the linear behavior 


C~-(7A-1) . 

TT 
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